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Abstract  Kerr  self-focusing  of  high-power  ultrashort 
laser  pulses  in  atmosphere  may  result  in  a  structure  or 
structures  of  high  intensity  that  can  propagate  over  long 
distances  with  little  divergence.  Filamentation  has  garnered 
significant  interest  in  the  nonlinear  optics  community  due 
to  its  unique  properties.  Salient  features  of  filaments 
include  a  central  region  of  intense  laser  power  (greater  than 
the  ionization  threshold  of  the  propagation  medium)  and  a 
low  temperature  plasma  column  that  lasts  up  to  nanosec¬ 
onds  in  duration  after  the  passage  of  the  laser  pulse.  Steel 
and  titanium  samples  are  ablated  by  filaments  and  by 
sharply  focused  sub-picosecond  laser  pulses.  We  then 
performed  metrology  on  the  samples  to  compare  the 
ablation  features  in  addition  to  modeling  of  the  plasma 
ablation  process.  Ablation  with  filaments  leads  to  a  wider 
range  of  material  responses  as  compared  to  ablation  with 
sharply  focused  pulse.  This  results  in  potential  complica¬ 
tions  for  applications  of  filament  ablation  that  depends  on 
the  rate  of  material  removal  and  spectroscopic  analysis. 

1  Introduction 

With  the  advent  of  nonlinear  optics  and  the  method  of 
chirped  pulse  amplification,  we  are  able  to  expose  matter  to 
extremely  intense  bursts  of  laser  light.  Of  particular  note  in 
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the  field  of  nonlinear  optics  is  filament  propagation  in 
atmosphere,  spawning  a  wide  array  of  studies  (see  review 
articles  in  11-3]).  The  filamentation  process  takes  advan¬ 
tage  of  air  molecules  that  would  disrupt  laser  propagation 
linearly,  however,  in  the  case  of  an  intense  sub-picosecond 
laser  pulse  can  act  in  a  nonlinear  manner.  For  example,  a 
1  mJ  collimated  laser  pulse  with  a  pulse  duration  of  100  fs 
(10  GW)  can  experience  Kerr  self-focusing  that  will  create 
an  increase  in  intensity  beyond  the  ionization  threshold  of 
air  (greater  than  lO^^W/cm^).  The  resultant  plasma  core 
defocuses  the  laser  light,  and  under  the  proper  conditions, 
the  competing  processes  create  a  filament  structure  that  can 
exist  for  hundreds  of  meters  14]. 

The  filament  structure  consists  of  a  plasma  column  of 
approximately  100  pm  in  diameter  15,  6]  surrounded  by  an 
energy  reservoir  as  shown  in  Fig.  1.  Of  the  initial  laser 
energy,  at  least  half  is  needed  in  the  reservoir  to  sustain  the 
intense  core  that  results  in  a  plasma  column  17].  The  energy 
reservoir  exists  on  the  same  timescale  as  the  initial  laser 
pulse  (sub-picosecond);  however,  the  plasma  column  takes 
much  longer  to  dissipate  the  energy  stored  in  both  the  free 
electrons  and  ions.  Thermalization  of  the  approximately 
1  eV  electrons  happens  well  after  the  laser  pulse  has  pas¬ 
sed,  and  the  recombination  time  is  on  the  order  of  nano¬ 
seconds  or  less  18]. 

Using  focused  ultrashort  pulse  lasers  (USPL)  for  abla¬ 
tion  is  attractive  due  to  the  precise  manner  in  which 
material  is  removed.  A  recent  review  of  USPL  ablation  as 
applied  to  micromachining  is  given  by  Cheng  et  al.  19]. 
The  hallmark  of  this  process  is  a  clean  crater  with  minimal 
spall  on  the  target,  whereas  nanosecond  (and  longer)  pulses 
will  leave  a  ragged  crater  and  material  strewn  along  the 
surface  110].  This  arises  from  rapid  transition  of  solid  target 
material  to  a  liquid  phase  via  superheating  during  the  laser 
pulse  111,  12].  Work  done  by  Kiselev  et  al.  113] 
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Fig.  1  Sketch  of  laser  propagation  in  which  Kerr  self-focusing  and 
competitive  plasma  defocusing  leads  to  the  production  of  a  filament 


demonstrates  the  inherent  benefits  of  the  concentrated  ul¬ 
trashort  energy  in  the  filament  structure  to  ablate  and  cut 
through  solid  matter,  termed  Filament-Induced  Laser 
Machining  (FILM).  As  opposed  to  complex  three-dimen¬ 
sional  machining  by  mechanically  moving  the  focusing 
lens  or  target,  a  filament  is  able  to  ablate  and  cut  along  a 
much  larger  depth  versus  the  Rayleigh  range.  The  micro¬ 
scopic  properties  of  filament  ablation  such  as  diameter  and 
volume  removed  show  less  than  30  and  50  %  relative 
standard  deviation,  respectively,  in  the  absence  of  geo¬ 
metric  focusing  over  distances  up  to  50  m  [14].  Related  to 
FILM  is  the  concept  of  filament-induced  breakdown 
spectroscopy  (FIBS),  essentially  utilizing  a  filament  to 
interrogate  the  chemical  make  up  of  a  target  or  a  coating  on 
a  target  in  the  manner  of  laser-induced  breakdown  spec¬ 
troscopy  [15,  16]  (LIBS)  but  at  remote  distances  [17-19]. 
In  our  study  of  how  the  focal  and  nonlinear  optical  con¬ 
ditions  impact  the  ablation  process,  we  have  defined  three 
regimes:  Sharply  Focused  Femtosecond  Laser  Machining 
(SFFLM),  Assisted  Filament-Induced  Laser  Machining 
(AFILM),  and  Lensless  Filament-Induced  Laser  Machining 
(LFILM).  SFFLM  is  achieved  using  what  we  consider  a 
short  focal  length  lens  that  reduces  the  possible  filamen- 
tation  regime  to  nearly  the  Rayleigh  range.  AFILM  uses  a 
long  focal  length  lens  in  addition  to  self-focusing  to  gen¬ 
erate  a  filament  ahead  of  the  geometric  focus;  the  filament 
extends  for  much  longer  than  the  Rayleigh  range.  Finally, 
LFILM  does  not  use  a  lens  but  relies  on  self-focusing  to 
generate  a  filament. 


2  Experimental  Setup 

For  the  experiment,  all  three  regimes  were  achieved  using 
the  same  USPL  system  but  with  different  focusing  princi¬ 
ples  as  shown  in  Fig.  2.  Our  laser  produces  pulse  energies 
of  up  to  25  mJ  in  100  fs  pulsewidth  at  800  nm  center 
wavelength  with  an  initial  beam  size  of  approximately 
1  cm  in  diameter.  The  SFFLM  regime  was  created  by  using 
a/^  =  25cm  lens  (similar  to  conventional  ablation  studies) 
and  the  AFILM  regime  using  a//  =  2  m  lens.  The  targets 
were  placed  before,  at,  and  after  the  focal  point  for  both 
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Fig.  2  Experimental  setup  showing  the  three  regimes:  a  Fast  focus, 
SFFLM;  b  Slow  focus  forming  a  filament,  AFILM;  and  c  collimated 
(no  optical  lens)  forming  a  filament,  LFILM 

SFFLM  and  AFILM;  the  positions  before  and  after  the 
focal  point  have  approximately  the  same  intensity  given  a 
linear  Gaussian  focus.  The  primary  method  we  used  for 
controlling  the  self-focusing  of  the  collimated  beam  in 
LFILM  was  by  changing  the  frequency  chirp  from  the  laser 
compressor.  As  such,  we  were  able  to  achieve  strong  self- 
focusing  in  the  LFILM  regime  as  close  as  3  m  from  the  end 
of  the  compressor. 

A  distinguishing  aspect  of  the  filament  structure  is  the 
long  region  of  intense  power  (necessarily  greater  than  the 
ionization  threshold  of  the  medium)  that  can  exist  for  many 
meters,  far  in  excess  of  the  Rayleigh  range.  To  examine 
how  that  is  manifested  in  the  ablation  phenomenon,  the 
target  was  placed  before,  at,  and  after  the  focus  position  for 
the  SFFLM  and  AFILM  cases.  Those  distances,  d  as  shown 
in  Fig.  2,  were  measured  from  the  lens.  In  LFILM,  the 
target  was  placed  at  positions  before,  immediately  after, 
and  well  after  the  position  where  self-focusing  is  pro¬ 
nounced.  Those  distances,  d,  were  measured  relative  to  the 
end  of  the  laser  system. 

For  each  ablation  point,  the  target  materials  were 
exposed  to  1,  10,  or  25  successive  pulses  (1  to  2  Hz  rep¬ 
etition  rate,  approximately  10  mJ  giving  100  GW  peak 
power)  to  examine  cumulative  effects.  The  materials  cho¬ 
sen  were  all  visually  opaque,  without  any  special  coatings, 
and  are  listed  in  Table  1  along  with  selected  material 
properties.  The  sample  surfaces  were  prepared  by  grinding 
to  a  finish  of  approximately  3  pm  peak-to-valley  for  the 
steels  and  6  pm  for  bulk  titanium.  The  samples  were  then 
examined  using  a  white-light  interferometer  (Newview 
6300,  Zygo  Corp.)  to  obtain  surface  profile  data.  The  root- 
mean- square  (rms)  variation  in  the  sample  surface  height  as 
measured  by  the  interferometer  was  0.06  pm  for  the  steel 
samples  and  0.30pm  for  bulk  titanium.  In  order  to  more 
closely  examine  the  shot-to-shot  variations  in  ablation,  we 
utilized  a  smooth  titanium  sample.  These  samples  were 
electron  beam  vapor  deposited  with  a  thickness  of  1  pm  on 
sapphire  substrate  to  a  rms  roughness  of  approximately 
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Table  1  Selected  properties  of  target  materials  at  298  K  [20-23] 


Material 

Main  composition 

Thermal 
conductivity 
(W/cm  K) 

Absorption 
coefficient  at 
800  nm  (|im^') 

Low  carbon 
steel  1018 

98.5-99.1  %  Fe 

0.803 

57.2 

Stainless 
steel  304 

62.5-74  %  Fe, 
18-20  %  Cr, 

8-12  %  Ni 

0.163 

47.0 

Titanium 
(bulk  and 
smooth) 

>99.0  %  Ti 

0.22 

52.1 

0.01  jim  and  placed  at  a  distance  of  11.8  m  from  the  laser 
for  LFILM. 


3  Results 

The  ablation  results  show  a  sharp  contrast  between  the 
short  focus  regime  (SFFLM)  and  the  two  filamented 
regimes  (AFILM  and  LFILM)  in  the  ability  to  remove 
material.  Figure  3  displays  the  results  from  the  low  carbon 
steel  target  and  is  qualitatively  representative  of  the  results 
seen  in  stainless  steel  and  bulk  titanium  (not  shown).  In 
Fig.  3a,  the  SFFLM  process  leads  to  a  sharp  crater  when 
the  target  is  placed  at  the  focal  point  of  the  25  cm  lens.  At 
the  SFFLM  focal  point,  the  Gaussian  intensity,  /,  approa¬ 
ches  2  X  lO^^W/cm^.  This  creates  a  well-defined  crater 


approximately  200  pm  wide  (about  5  times  the  spot  size) 
and  a  lip  with  a  peak-to-valley  height  difference  of  about 
5  pm  with  approximately  5  x  10“^mm^  (0.39  |ag)  of 
material  removed.  In  this  sharply  focusing  regime,  the  laser 
intensity  surpasses  the  ionization  threshold  of  air,  on  the 
order  of  10^^  —  lO^^W/cm^,  leading  to  breakdown;  this 
region  extends  3  cm  on  either  side  of  the  focal  point. 
However,  we  can  see  from  the  average  cross  sections  (over 
1 1  lines  of  data)  in  Fig.  4a  that  placing  the  target  6  cm  on 
either  side  of  the  focal  point  {d  =  19cm  or  31  cm,  spot  size 
of  2.1  mm,  /  =  6.4  X  lO^^W/cm^)  results  in  little  ablation. 
Similar  results  are  seen  in  stainless  steel  (Fig.  4b)  and  bulk 
titanium  (Fig.  4c). 

Switching  to  the  AFILM  regime  with  a  2-m  focal  length 
lens,  it  produces  a  marked  change  in  the  ablation  results. 
As  shown  in  Figs.  3b  and  5a,  a  less-defined  crater  is  seen 
on  low  carbon  steel  at  the  focal  point  though  the  affected 
area  is  almost  twice  as  wide.  The  depth  in  the  center  of 
ablation  area  is  shallower  than  in  SFFLM  though  this 
results  in  roughly  the  same  volume  of  material  removed. 
This  is  despite  the  peak  intensity  being  almost  two  orders 
of  magnitude  less.  Similarly,  the  stainless  steel  target  in 
Figs.  5b  and  6  shows  a  clear  deviation  from  a  circular 
feature  in  the  ablated  area  at  2  and  2.47  m  from  the  lens 
though  the  spot  size  and  intensity  are  similar  to  the  19  and 
31  cm  positions  in  SFFLM  (Fig.  4a),  demonstrating  that 
we  are  in  a  nonlinear  optical  regime.  This  structure  may  be 
the  result  of  the  laser  pulse  entering  the  multiple  filament 
regime  where  individual  filaments  can  form  from  the  initial 


Fig.  3  False  color  height  data 
of  low  carbon  steel  target  after 
25  pulses  with  a  SFFLM  A  = 

25  cm  at  =  25  cm,  b  AFILM// 
=  2  m  at  J  =  2  m,  and  c  LFILM 
at  J  =  10.4  m 


Fig.  4  Average  cross  section  of  low  carbon  steel  (a),  stainless  steel  (b),  and  bulk  titanium  (c)  targets  along  the  direction  of  laser  polarization 
after  25  pulses,  with  SFFLM  A  =  25  cm.  Colored  lines  denote  different  target  distances,  d,  relative  to  the  lens 
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Fig.  5  Average  cross  section  of  low  carbon  steel  (a,  d),  stainless  steel  (d-f)  LFILM,  no  lens.  Colored  lines  denote  different  target  distances, 

(b,  e),  and  bulk  titanium  (c,  f)  targets  along  the  direction  of  laser  d,  relative  to  the  lens  (a-c)  or  laser  system  (d-f) 

polarization  after  25  pulses,  with  (a-c)  AFILM  //  =  2  m;  and 


Fig.  6  False  color  height  data 
for  stainless  steel  target  after  25 
pulses  with  AFILM  //  =  2  m  at 
a  J  =  1.5  m,  b  J  =  2  m,  and 
c  d  =  2.47  m  from  the  lens. 
Scales  are  identical  for  each 
image 


2.5  pm 


0 


Fig.  7  AFILM  comparison  at  focal  point  (J  =  2  m)  of  a  single  pulse 
burn  paper  {white  2  mm  grid  lines),  and  b  25  pulses  on  low  carbon 
steel 

pulse  [24-26].  At  the  peak  power  of  our  pulse  (100  GW), 
we  could  expect  up  to  10  individual  filaments  to  form  [27]. 
However,  we  observe  on  bum  paper  (Fig.  7a)  a  non-cir¬ 
cular  area  that  may  result  from  multiple  filaments  that 
would  correspond  to  the  complex  ablation  pattern  seen  in 
Fig.  7b.  It  is  important  to  note  the  sharp  topographic  var¬ 
iation  between  Fig.  6a  and  c  despite  being  almost 


equidistant  from  the  focus  demonstrating  a  substantial 
change  to  the  profile  of  the  laser  pulse.  However,  the  peak- 
to-valley  surface  height  varies  little  at  the  different  values 
of  d,  in  comparison  with  that  seen  in  SFFLM  (see  Fig.  5a- 
c).  The  bulk  titanium  target  showed  less  susceptibility  to 
material  removal  at  any  distance  under  AFILM,  but  rather 
showed  a  pronounced  visual  change  to  the  surface  (i.e., 
darkening).  This  result  is  in  sharp  contrast  to  the  observed 
SFFLM  ablation  material  removal  and  crater  depth  in 
Fig.  4c  when  compared  to  the  steel  targets. 

Removing  the  lens  and  adjusting  the  frequency  chirp, 
we  then  created  a  filament  based  on  self-focusing  without 
geometric  focusing.  In  this  LFILM  regime,  we  start  to  see 
strong  self-focusing  at  3  m  from  the  end  of  the  laser  system 
and  strong  continuum  generation  by  5  m.  Figure  3c  shows 
an  ablation  area  on  low  carbon  steel  that  is  not  as  sharply 
defined  as  in  the  SFFLM  (Fig.  3a)  and  AFILM  (Fig.  3b) 
regimes  even  at  10.4  m  from  the  laser.  The  ablation  area  is 
even  larger  than  SFFLM  and  AFILM,  possibly  due  to 
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higher  pointing  instability  from  a  lack  of  geometrical 
focusing,  with  the  ablation  depth  less  and  only  about  half 
the  volume  of  material  removed.  The  remarkable  result  of 
LFILM  is  the  reproducibility  of  the  ablation  features  of  the 
steel  targets  over  a  wide  variation  in  d  as  shown  in  Fig.  5d, 
e.  As  in  the  AFILM  case  for  the  steels,  the  bulk  titanium 
target  seen  in  Fig.  8  was  much  less  affected  by  LFILM  in 
terms  of  ablation  and  material  removal,  however,  showed  a 
significant  darkening  of  the  surface  [28,  29]. 

We  examined  the  smooth  titanium  targets  in  LFILM  and 
see  in  Fig.  9a  a  crater  size  of  approximately  0.5  mm  in 
diameter  in  agreement  with  the  bulk  substrate  target 
(Fig.  8).  A  central  peak  of  material  rises  above  the  ablation 
crater  floor  in  the  single  ablation  spot  in  Fig  9a,  and  the 
maximum  transverse  beam  movement  can  be  measured  by 
quantifying  the  distance  between  peaks  produced  by  two  or 
more  sequential  laser  pulses.  In  this  example,  the  peak-to- 
peak  distance  is  less  than  0.4  mm  (Fig  9b).  We  expect  the 
most  variation  in  transverse  beam  movement  in  LFILM 
due  to  a  lack  of  geometrical  focusing.  As  can  be  seen  in 
Fig  9b,  the  movement  of  the  laser  pulse  can  have  both 
parallel  and  perpendicular  motion  relative  to  the  laser 
polarization  (in  this  figure  vertical).  Increasing  the  number 
of  pulses  results  in  a  preferential  elliptical  region  of  abla¬ 
tion  of  approximately  0.78  mm  diameter  parallel  to  the 
laser  polarization  direction  and  0.63  mm  diameter 


Fig.  8  LFILM  of  bulk  titanium  target  at  J  =  10.4  m  after  25  pulses: 
a  raw  monochrome  image  and  b  false  color  height  data 


Fig.  9  Smooth  titanium  target  after  a  1  pulse  and  b  2  pulses  LFILM 
at  J  =  11.8  m  from  the  laser  system,  false  color  height  data 


perpendicular  after  100  pulses  (data  not  shown).  Examin¬ 
ing  the  detectable  peak-to-peak  movement,  the  average 
over  all  the  samples  was  0.21  mm  with  a  standard  devia¬ 
tion  of  0.13  mm. 


4  Plasma  ablation  simulations 

In  order  to  understand  the  ablative  properties  of  the  plasma 
component  of  the  filament  in  FILM,  a  plasma  ablation 
model  was  developed  to  represent  the  possible  effect  the 
lingering  plasma  column  may  have  on  the  target  surface. 
The  plasma  ablation  model  used  here  has  successfully 
modeled  atmospheric  ablation  of  pulsed  plasma  thrusters, 
plasma-propellant  interactions,  and  arc  discharges  and  is 
fully  detailed  in  [30].  The  model  is  based  on  the  concept  of 
coupling  Anisimov’s  vacuum  laser  ablation  theory  to  a 
hydrodynamic,  non-equilibrium  layer  to  relate  the  target 
properties  to  the  bulk  plasma  properties  [31].  In  order  to 
close  the  model’s  system  of  equations,  the  vapor  pressure 
of  the  target  is  required.  A  Clausius-Clapeyron  relation 
was  fitted  using  data  from  room  temperature  [32],  having 
found  that  the  final  temperature  of  the  target  does  not 
change  significantly.  Finally,  the  model  consists  of  solving 
a  single  algebraic  equation  numerically  to  determine  the 
ablation  rate  of  the  target. 

The  ablation  rate  is  dependent  on  the  instantaneous 
surface  temperature  of  the  target  and  may  change  with 
time.  Thus,  the  ablation  model  (detailed  description  given 
in  Chapter  4  of  Ref.  [30])  was  coupled  to  an  axisymmetric 
two-dimensional  thermal  conduction  model  to  represent 
the  target  depth  (see  Fig.  10).  A  finite  plasma  column 
radius  of  50  pm  was  assumed,  while  the  computational 
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r 
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Fig.  10  Simulation  layout  used  with  filament  impinging  on  metal 
target  in  axisymmetric  model 


Table  2  List  of  material  properties  at  298  K  and  1  atm  used  in 
simulations  [20,  21] 


Material 

Density 

Specific  Heat  at  Constant 

Heat  of 

g/cm“^ 

Pressure  J/kg  K 

Vaporization 

J/kg 

Iron 

7.86 

450 

6.09  X  10® 

Titanium 

4.51 

523 

8.88  X  10® 
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domain  had  a  radial  extent  of  five  times  the  column  radius. 
Target  depth  was  half  the  radial  extent  making  it  func¬ 
tionally  semi-infinite  for  the  timescales  considered  here. 
The  thermal  model  requires  thermodynamic  properties  (see 
Table  2)  for  the  target  such  as  thermal  diffusivity,  heat  of 
vaporization,  and  density.  For  temperature  dependent  val¬ 
ues,  room  temperature  data  was  used. 

Some  assumptions  on  the  state  of  the  plasma  column  had 
to  be  made  due  to  lack  of  experimental  data  probing  the 
filament.  Total  heat  flux  reaching  the  target,  bulk  plasma 
number  density,  and  bulk  plasma  temperature  are  required 
inputs  to  the  model.  For  these  parameters,  values  from 
typical  atmospheric  density  plasma  ablation  were  used 
and  were  selected  to  yield  the  greatest  calculated  target 
ablation.  These  parameters  are  lO^W/cm^,  lO^^crn”^, 
and  1.5  eV,  respectively.  Targets  simulated  were  titanium 
and  iron  (similar  in  composition  to  the  low  carbon  steel 
target  material). 

Based  on  known  properties  of  filaments,  it  was  assumed 
that  the  plasma  column  will  only  persist  in  the  atmosphere 
for  1  ns  after  the  laser  pulse.  Unsurprisingly,  the  model 
shows  very  little  temperature  rise  in  the  target  and  no  target 
ablation  for  either  titanium  or  iron.  Yet  the  model  still  tells 
us  something  important  about  filaments,  namely  that  the 
plasma  column  alone  is  not  responsible  for  any  of  the  fil¬ 
ament-exclusive  shapes  discussed  in  the  previous  sections. 
Even  by  increasing  the  plasma  duration  to  100  ns,  the 
target  temperature  change  does  not  result  in  ablation.  Thus, 
the  effects  of  the  ablation  on  the  targets  shown  above  are 
due  to  interaction  with  the  initial  laser  pulse,  not  post-pulse 
interaction  with  the  plasma  column.  The  plasma  column 
may  have  some  affect  on  the  filament-target  interaction, 
but  the  modeling  results  discussed  here  indicate  that  that 
effect  is  not  thermal  or  ablative  in  nature. 


5  Conclusion 

The  filament’s  unique  behavior  to  create  an  intense  and 
narrowly  diverging  pulse  of  power  with  a  plasma  column 
over  long  distances  remains  an  intriguing  method  to  ablate 
solid  materials.  In  contrast  to  using  a  short  focal  length 
lens,  a  filament  can  ablate  and  mark  materials  over  a  much 
longer  distance  (10-100  m).  As  evidenced  by  removing  the 
focusing  lens,  we  see  nearly  reproducible  ablation  char¬ 
acteristics  over  many  meters  though  these  are  not  as  well 
defined  as  when  a  short  focal  length  lens  is  used.  The 
drawback  of  filament  ablation  is  a  smaller  amount  of 
material  removed  that  could  be  of  importance  to  an 
application  like  FIBS  with  the  additional  variables  of 
focusing  paradigm  and  substrate  material.  Using  a  short 
focal  length  lens  is  clearly  superior  if  the  goal  is  to  perform 
a  precise,  well-defined  cut  with  the  highest  rate  of  material 


removal.  Ablation  with  a  filament  shows  promise  when  the 
target  is  at  varying  or  far  distances  from  the  laser.  How¬ 
ever,  the  results  shown  here  lead  to  further  investigation 
into  understanding  how  target  material  properties  lead  to 
substantial  differences  in  results  from  filament  ablation. 
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